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ABSTRACT 
The objective to design a PID Controller with negligible 
parameter interdependency was carried out. As a result the optimum 
tuning graphs developed by A.M. Lopez, based on integral of error 
criteria, ideal PID algorithm and steepest descent computing technique 
can be used in conjunction with this near ideal PID controller. Reset 
inhibiting and other convenient features are also incorporated to make 


this controller more versatile. 
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Chapter I 
Introduction 

The Conventional Continuous PID Controller 

In its early days, the design and application of the 
conventional continuous PID controller was somewhat of an art [1]. 
The choice of the PID combination itself was guided by intuition 
and experience. In those early days, the problem caused by dead 
time delay and physical limitations were great obstacles to 
analytical study of the system [2]. Hence, it is not surprising 
to find that even the tuning methods were derived by trial and 
error. Quite a number of tuning methods are in use. To name a few: 
Ziegier-Nichols method, Harriot's method, 3-C method etc. Most 
of these methods were based on the quarter decay ratio (or some other 
arpLerary constraints which fack»ssolid mathematical justification) 
as the performance criterion. 

To complicate the matter further, most of the conventional 
PID controllers have their three parameters interrelated [3]. Thus, 
even if the correct parameter settings were known, it could well 
take hours to settle down to the final adjustments. It is therefore 
obvious that optimal tuning of a controller is quite a difficult 
problem (By optimal tuning is meant that tuning which gives the 
best performance, judged by a certain mathematically sound per- 


formance criterion). 
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Possible Improvements 


The optimal tuning problem should not be avoided, because 
it can be shown that efforts devoted to optimal tuning can bring 
spectacular results [3]. 

The tuning difficulties outlined above actually initiated 
the development of the digital version of PID controller, which 
has independent parameters. It was shown in reference 4 that as 
sampling frequency increases the overall performance of a digital 
PID controller improves. In the limit it approaches the performance 
of a continuous controller. This means continuous PID controller 
would be superior but for its difficulties in optimal tuning. 

itjwas alsoeproved by M, Athans and K.1. Parker [5, 11]% 
that proper Pee of cost functionals for the general tracking 
problems, and application of optimal linear regulator theory to the 
optimal tuning problem indeed yield a PID type of controller as the 
solution. 

The above arguments with respect to continuous systems and 
optimality provide sufficient grounds to warrant taking a second 
look at continuous PID controllers. 

A.M. Lopez in his Ph.D. dissertation [6] had indicated the 
shortcomings of the quarter decay ratio as a performance criterion, 
especially if applied to higher order plants. He also indicated the 
desirability of using the integral of error as a performance index. 


He used a steepest descent computer program to search for optimal 
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settings for plants with dead time delay (using the integral of 
error as performance criterion and assuming an ideal algorithm with 
no lag and independent parameters). 

A question arises as to whether the ideal algorithm can be 
closely approximated. (In real controllers, no lag derivative 
action can not be achieved because of noise). If a close approx- 
imation can be achieved, one can use the results provided by Lopez 
as the optimal settings for a system with integral of error as the 
performance criterion. Assuredly this three fold improvement 
(better performance criterion, independent parameters and computer 
derived results) will bring about spectacular results when compared 
to the previous trial and error situation. 

In both pneumatic and electronic controllers with purely 
passive networks, the parameter interrelating situation is unavoid- 
able [3]. This problem can be solved by using active networks. 
With the advent of integrated circuits, independence of parameters 
can be easily achieved. Also, due to the negligible D.C. drift of 
integrated circuits, one can avoid complicated D.C. to A.C. and A.C. 
to D.C. conversions as was necessary in conventional electronic 
COULTOMLEES. [7 5.0 | AS CeLtaln amountsoLnelagemust exist in the 
derivative action. Otherwise the gain at high frequencies will be 
intolerably large. However, the objective in this thesis project was 
not to build a completely ideal controller, but rather to develop 
as close an approximation to it as possible. A.M. Lopez's ideal PID 
algorithm will be shown in Chapter IV. PID stands for proportional, 


integral and derivative. 
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The Objective 


The primary objective in this thesis, therefore, is to design 
a continuous PID controller with an algorithm closely approximating 
the ideal one cited by Lopez, such that all the tuning graphs 
developed by Lopez will be directly applicable. 

Other features that will be incorporated in the controller 
are: integral inhibiting action, direct and indirect action, 
manual control, bumpless transfer, digital and analog set point 
adjustment. All these are included to make the controller more 


versatile and suitable for industrial usage. 
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Chapter II 


Theoretical Background 


The conventional PID controller is based on the simple idea 
of negative feedback. It is well known that the linear regulator 
theory requires only state variable feedback [9]. In essence, the 
plant output is fed back and compared with the set point (input). 
The difference between them, the error signal, will be manipulated 
in a certain fashion to generate the control action. In turn the 
control action will force the plant output to move in a direction 
such as to reduce and eventually eliminate the error signal. There- 
fore, the sole purpose of the controller is to provide the desired 
manipulation of the error signal. The position of a controller in 


a control loop-is shown in figure 2-1. 


CONTROLLER 








CON TROL 
OUTPUT OUTPUT 





Figure 2-1 The position of the controller ina 
control loop. 
Proportional Action (P action) 
The simplest way to manipulate the error signal is to multiply 


it by a gain Ko as the following: 

















Sp Oios 
Ae), (Pooh sas ac Sane el reblog. 10> 017 Lenal tumeaiia a mate , a 





; t ; / af 
PHEBAT WHSOEL ac eile wath, bisy=e2 MTG Le Dp ‘pul 


7 _ 
_ 
f 


7 





— 


sage <4hse aoe al? side, BSRaaee ne na. fet re eet 


a 


P ee) al be ivf ,ai hae Set er ee 
af Tye, a L258 ie I Je a eid a YSIS a3, Magee aso ee ' 


; isiD & at sven o3.9uqtee oa ly Sao Saran Lie. 3 ee 


- ru J 
HM lamgte 25215 31 ibis yl ieogtsey- bigoasulsT-se ad doe 
<< - - — 


~ 


se am pee terions ods teins 


a2, z - e 





m(t) = K e(t) 
| (2-1) 


M(S) 


il 


K ECS) 


where m(t) is the control action 
e(t) is’ the error signal 
RK. is a gain constant 
and § is the Laplace operator. (assuming all initial conditions are Zero) 


The proportional action is shown in figure 2-2. 





Figure 2-2 The proportional action in a control loop. 
From Figure 2-2 it is obvious that the control action m(t) is directly 
proportional to gain K. and error signal e(t). A higher value of 
m(t) will reduce the error faster. However, for a "proportional only" 
controller there will always be some steady state error (offset), 
otherwise the control action RCD would be zero. Although higher 
gain will produce a smaller offset, it might also lead to a reduced 
degree of stability in the control loop. 
Proportional Plus Integral Action (PI action) 

The integral action is primarily introduced to cancel offset. 


By definition it is the action which integrates the error signal. As 
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error persists the integrated value will grow with time, and 
eventually will be large enough to eliminate offset completely. 


The formulation, o& PL action is: 


K 
m(t) = K_ e(t) ae e(t) dt 
i 
K E(S) 
M(S) = K E(S) + - : 
aL 
1 
=K (1+ Ts FC) (222) 


where T, is the integral time constant. Or in block diagram form in 
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Figure 2-3 PI action in a control loop. 
Note that if e(t) equals a constant and t=T, in equation (2-2), the 
integral part will be equal to the proportional part. Therefore, the 
meen time constant qT, is the time needed to produce an integral 
action equal to the proportional action. This situation is depicted 


in figure 2-4. 
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Figure 2-4 The relationship between integral time constant 

qT, and control action. 
At the very instant when disturbances or load changes of a process 
occur, the integral action is very small, because the time elapsed 
Usestitleveryeshort.) At this early time the major action is the 
proportional action. Later the integral action takes over and forces 
the error to diminish and eventually vanish as time goes on. As the 
error approaches zero so does the proportional action. However, if 
the values of 1/T, is too large the integrator will charge up too 
much and cause the output response to overshoot the set point. Only 
after overshoot has actually occurred (the error signal will change 
sign) will the integrator start to discharge. Actually, too much 
integral action will decrease the stability of the loop. 


The frequency response of PI action is shown in figure 2-5. 
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Figure 2-5 The frequency response plot of PI action. 

Note that at higher frequencies proportional action predominates. 
This is in accord with the previous analysis; because the initial 
transient of the response consists mainly of higher frequency 
components. On the other hand, as steady state is approached the 
output contains primarily very low frequency components and hence 
mainly the integral action is of importance. 
Proportional Plus Derivative Action (PD action) 

Derivative action is needed whenever substantial time lags 
exist in, the: control loop... This actionsis proportional to the) rate 
of change of error signal, and hence is predictive in nature. The 


PD action can be described as: 


m(t) = K.e(t) + KT, ic) 
M(S) = K,E(S) + KT, 8 E(S) 


Ka aE Ty S) E(S) (2-3) 
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where T is the derivative time constant. Its block diagram form 


is shown in figure 2-6. 


Ke (+ 1d 5) 
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Figure 2-6 PD action in a control loop. 
At the very instant when a disturbance occurs, the 
derivative action is in the same direction as the proportional action 
and therefore fulfills its anticipatory function. After the initial 
period the derivative action will become opposite in sign with 
respect to the proportional action. This, in effect reduces the control 
output, slows down the process and therefore increases the stability 


of the control loop. This situation is shown in figure 2-7. 
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sudden change of set point. 

Usually, when derivative action exists, the proportional 
gain K. can be increased because of the increased stability. Under 
steady state the derivative action is zero . The above description 
can also be derived from the frequency response plot of PD action 


as shown in figure 2-8. 
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Figure 2-8 The frequency response plot of PD action 
If e(t) = Ast, a ramp function (A is a constant), the control 


action will ber: 


7d e(t) 


m(t) = Kie(t) + KT Ar 


=I Ree Att ke Tee A (2-4) 
Cc Cc 


The derivative part of control action is equal to KTgtA- When 


m(t 4) = eel ant KTyA (2-5) 


d 


That is to say, when t = T the proportional part will be equal to 


(al? 


the derivative part as shown graphically in figure 2-9. 
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Figure 2-9 The relationship between derivative time 
constant qT, and PD action 
It should be clear from figure 2-9 that the derivative action is Ta 


units of time ahead of the proportional action. 


Proportional Plus Integral Plus Derivative Action (PID action) 





The formulation of the three mode controller is: 


K 
avers (Fe 2S 
m(t) Ke eCl) = i, e(ryedt+ KT 


de(t) 
dade 


U 


NUE): SPS GE GS ep | 3 T 4S) E(S) (2-6) 


i 
ae 
ae 


The block diagram representation of equation (2-6) is shown in figure 2-10. 





Figure 2-10 Block diagram of a PID (three mode) controller 
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If e(t) is again a ramp function the control action will be as shown 


inetd Sure w 2 = le 


m (+) 





Figure 2-11 Three mode controller action with ramp input 


The frequency response of the three mode controller is shown in 


figure 2-12. Ww, and Wz in Figure 2-12 are derived from eqaution 2-6. 
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Figure 2-12 Frequency response of a PID controller 
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In figure 2-12 the high frequency end exhibits a positive slope. 

This positive slope represents the effective derivative action and is 
intended to take care of errors with primarily higher frequency 
components. Similarly the effective integral action is represented 
by negative slope at the low frequency end and will deal with steady 
state errors which consist essentially of low frequency components. 


The manipulation of T, and qT, will shift the two corner frequencies 


d 
of the response horizontally, while a change in K. value will only 
shift the whole curve up and down. 

Generally the integral action will decrease stability, reduce 
offset to zero and will increase the period of any possible ringing. 
The derivative action will increase stability and decrease the period 
of ringing. Therefore, when derivative action exists one can use 
higher values of K. and 1/T,. Unfortunately, the derivative action 
is also very sensitive to noise. 

References 1, 2 and 10 give detailed discussion on this subject. 
Modern Point of View 

In the past, the application of optimal linear regulator 
theory to industrial control problems led to PD type of linear feedback 
systems. As recently as 1972 M. Athans [5] showed that, in choosing 
cost functionals, consideration has to be given to three factors. 
These factors are: the cancelling of steady state error, the boundedness 
of control output and the independence of control output from nominal 
plant parameters. The resultant control algorithm for a first order 
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in the same year, expanded Athan's method to cover the general problem 
of tracking m variables with a nth - order system. For example, the 


following control loop will result for a second order plant (figure 
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Figure 2-13 The optimal PID combination for the second order plant 
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action. Therefore, the works of Athans and Parker can serve as a 
solid a OCIA basis for three mode controlling which was 
implemented long ago by intuition and experience. 

In this thesis, however, it is not intended to design an 
optimal controller in the strictest sense because great effort is 
required merely to find a specific PID combination for a particular 
plant. Instead, to be more general, a fixed format PID controller 
with independently variable parameters will be developed. Such a 
controller, with easily adjustable parameters may then be applied 
to a wide variety of plants. 

Integral Inhibiting Action (Reset Inhibiting) [12,13] 

By using integral action, the offset inherent in proportional 
only control will be eliminated. However, when PI action is applied 
to processes requiring either start up, shut down, large changes in 
control point or major load changes, large overshoots are generally 
present during the transient, because of a condition known as integral 
wind-up. In critical processes such as batch-sintering, batch brazing, 
diffusion, crystal growing etc., overshoot may not be tolerated 
because it can affect the quality of the product. Therefore many 
critical processes using conventional controllers are started manually 
to avoid overshoot. 

At the beginning of a sudden change, the error signal is 
large. The integral action will integrate this large error signal. 


As time goes on it probably will saturate the integrator. Consequently 
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the controller output will be at its maximum or minimum values. 
This situation is called integral wind-up. The integral action will 
get out of saturation only when the error signal changes sign i.e. the 
controlled variable cross the set point (overshoot). Once overshoot 
occurs, the integrator will start to discharge from the value it 
stored up during the initial moments. It might take some time to 
discharge and eventually bring the controlled variable back towards 
the set point. 

One scheme devised to remedy this situation is to use 


derivative action. Consider the responses in figure 2-14: 


eGo) Kc e(%) 





Kepfete) dt KeTade ) 


Figure 2-14 The various responses to a set point change 
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Except in the initial instant, the derivative action is in opposite 
direction to both P and I actions. Therefore the D action can cancel 
part of the I action. As a result, the overshoot will be reduced to 
some extent. Nevertheless, this method is not adequate to solve the 
problem, especially when the steady state controlled variable is 

only a fraction of the maximum controller output. (i.e. the integrator 
has a great deal of range to charge up to a large value) 

The reset inhibiting scheme to be used in the proposed design 
is derived from the following reasoning: (1) Both integral and 
derivative action are aimed primarily at on-line control where integral 
action can eliminate offset and derivative action yields quick response 
to disturbances. (2) In a practical controller there will always be 
a limit to the range of output (output bounded). Initially the error 
is so large that proportional action alone would probably demand 
maximum or minimum controller output. Therefore, the reset action 
can actually be shut off initially. 

To decide the proper moment and amount at which to re-introduce 
integral action, the following study of figure 2-15 is required. In 
figure 2-15, the error signal at point A will cause the proportional 
action to be either 70% or 62% of the maximum output,depending on the 
gain setting Ke SO, Longe as verrotreexicts. Plaction: willebe ditrerent 
from its steady state value. The value of P action will depend on 
both the size of error and the P action gain setting. If the I action 


is to assist the P action, it should be in the same direction as 
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Figure 2-15 The proportional action versus error signal with 
different gain Ko as parameter 

P action. In the extreme, when P and I action combined give out 

the maximum or minimum value of output, the response will be the 
fastest. Therefore, it is desirable that the proportional and integral 
action combined have the capabtlity of producing output extremes. 
For example, if the error requires the P action to produce 75% output, 
the I action should be capable of contributing the remaining 25% of 

the output range in the same direction. During steady state the P 
action becomes zero. Hence the integral action alone should be capable 


of giving the output extremes. It should be realized that reset 
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wind-up is no longer of concern because the reset action is initially 
inhibited. 

To conclude: whenever the error signal is within the 
proportional band (and therefore the P action is within the output 
range) the I action will be allowed to attain such a value that the 
combined P and I action equal 100% of output range. This can be 
achieved by constantly monitoring the sum of P and I action. Whenever 
the sum exceeds the output limit, the integrator will be completely 
discharged, in other words, inhibited. With this arrangement, the 
integrator becomes self regulating. It adjusts its own operating 
range to precisely that amount needed to produce output extremes. 
Whenever the extremes are exceeded, it inhibits its own operation. 

Besides the elimination of the reset wind-up problem, other 
benefits derived from this scheme are: 

(1) there will be no start-up problam, and 

(2) the derivative action is free to carry out its original 


purpose. 
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Chapter III 

The Design of a Near Ideal PID Controller 

With the design objective and its theoretical background 
already laid out in Chapter I and II, the design work itself may now 
be discussed. 
System Description 

For any real device, physical limitations are unavoidable. 
Hence, an understanding of the operational limits of the controller 
and the control loop is needed. In figure 3-l, a particular control 


loop and its ,signal operating ranges are shown. 






PHYSICAL VARIATIONS 





Gh) 


c(t) (i-S5V) 





Figure 3-1 A particular control loop and its signal operating 


ranges 


The distance between a controller and the controlled plant could be 


quite large. In order to prevent the electrical signal from being 
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affected by the distance, the controller output port is designed to be 

a current source. Actually the controller internal signal of 1 to 5 

volt is converted into 10 to 50 mA current. This is because the majority 
of industrial transducers operate in the 10 to 50 mA range. 

The sensor within the plant will provide the voltage 
representation of the controlled variable. This voltage output eX Gad) 
will be subsequently transformed into 10 to 50 mA range of current 
The reason for the transformation is again the possible long distance 
between controller and plant. Upon arrival at the controller, the 
output c(t) (in current form) will be converted back into a ito 
5 volt range of voltage. The set point voltage is also in the range 
of 1 to 3) volte|tee5]- 

The Ba onceed PID controller will have the following ranges 


for its parameters: 


Proportional gain K. 0 SiR. ss Al) 

Integral time T, 0 << 10 repeats/min 
af 

Derivative time Ty 0 sho < 10 minutes 


As will be shown in Chapter IV, these parameter ranges are adequate 
for most industrial requirements. 
Power Supply 

Since integrated circuits will be used in the controller and 
a larger signal swing is desirable, + 15 volt is used as the power 


supply voltage. A + 40 volt supply is also needed for the output 
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current drive. The + 40 volt supply was selected for convenience in 
this proposed design. If necessary a 90 volt supply can be used to 
allow series connection of transducers. 

LG eeprecen ely quite common to rack-mount controllers. 
Hence, it is convenient to have one set of external power supplies to 
provide power for several controllers and their accessories. In this 
proposed design, only external power supplies will be used. The 
required + 15 volt and + 40 volt supplies are drawn from standard 
laboratory power supplies. The auxilliary + 6 volt and + 3 volt 
supplies for logic circuits are derived from the + 15 volt supplies 
by voltage regulators and voltage divider circuits. The complete 
circuit details for the power supplies are shown in figure 3-18. 
Manual and digital set point adjustments 

The manual set point voltage can be easily generated as shown 


in figure 3-2. 
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Figure 3-2 The generation of manual set point variations 
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The + 1 to + 5 voltage range is derived from the potentiometer 
and fed into the positive input of the operational amplifier. The 
Op-amp here ere as a buffer and a voltage source. 

The digital set point adjustment was designed into a part- 
icular format. This particular design is intended only to demonstrate 
FHespossipility of digital computer control. | Othergforms=otydigital 
set point adjustment are possible. 

For this particular design, the digital set point range of 
1 to 5 volt is divided into one hundred evenly spaced steps, such 
that each step equals 40 mV. [(5-1)volt/100 = 40 mV]. The command 
signal, asking for any quantized value between + 1 to + 5 volts, is 
coded in binary format. To represent 10055 in binary, one needs at 
least seven bits. For example, the binary 1100100 equals the decimal 
100 and will represent a + 5 volt output (+ 1 volt + 40 mV X 100)- 

The devices chosen to convert binary codes into voltages are 
a Fairchild's yA722 10 bit current source and an op-amp as shown in 
figure 3-3. The resistors marked as Ro in figure 3-3 are selected, 
such that each successive transistor in the pA722 device will deliver 
twice as much current as the previous one. By applying a binary 1 
to any logic gate in the »A722 device, the gate will be closed and 
the transistor below that gate will draw an integral number of current 
units from the following op-amp. As specified in the yA722 


specification sheets, the maximum current when all bits are on, 
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Figure 3-3 Using wA722 and an op-amp to convert a 
binary code into a voltage 
is 2560 uA. Taking the current through the least significant bit 
transistor as one unit, the maximum atazae when all bits are on, 
corresponds to 1023 units (N = y 2° = 1023). Therefore, each unit 
equals 2.5 pA (2560 EVRUEER REE ON Since only seven bits are 
needed to represent 100 steps in binary, the seven most significant 
bits in the A722 device were chosen. This way the ywA722 will be 


able to deliver a heavier output current. The least significant bit 


in the A722 device will now be the 7th bit, which delivers a current 
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8 times (23 = 8) the current of the 10th bit. Therefore, the 7th 
fransistorewill draw 20 wA (2.5 wA x 8), which is the new current unit. 
For each step to correspond to 40 mV, the 20 wA drawn through the op- 
amp should produce 40 mV. Consequently 2K2 (2K2 x 20 yA = 40 mV) 
was chosen to be the feedback resistor value for the op-amp. 

Manipulation of the binary digits is accomplished through 
mechanical switches on the front panel. But these mechanical switches 
are not directly connected to the gates of the uwA722. Instead the 
output of each mechanical switch goes to a 'D' type flip-flop. The 
'D' type flip-flop has to wait for a read pulse before it transfers 
its content to the gates of the pA722. The purpose of the above is 
to allow time for manipulation of the mechanical switches, without 
changing the previous set point level until the read-in switch is 
activated. The read-in switch is a normally open mechanical switch, 
which will give a + 3 volt output when pushed on. At the release of 
the read-in switch, the trailing edge of the + 3 volt pulse will 
trigger a monostable multivibrator, which in turn will produce a one 
shot pulse. The leading edge of the one shot pulse will act as the 
read pulse to the 'D' type flip-flop. 

The width of the one shot pulse is calculated to suit the 


particular 'D' type flip flop chosen. 
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Aepedestalsof | volteis added to the output to elevate the 
set point level ro the + 1 to + 5 volt range. The complete circuit 
is shown in the digital set point section of figure 3-18. 

Error Signal Generation 

By definition the error signal is the difference between the 
set point level and the controlled variable. Assuming that the 
feedback signal is in current form, then a conversion to a voltage 
signal is needed. After the current to voltage conversion the plant 
output and the set point voltage are subtracted from each other by 
feeding them into the opposite polarity inputs of a differential 
amplifier. Therefore the output of the differential amplifier is 
the error signal. The circuit to generate error signal is shown 


in figure 3-4. 
-(L) SET POINT 
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29 
The double pole double throw switch in figure 3-4 serves the 
purpose of generating either + e(t) or - e(t). This switch therefore 
enables one to select direct or indirect action. 


The Controller Proper 


The ideal PID algorithm in mathematical formulation is: 


—— = K (1+ 


M(S) _ i aw 
E(S) c Ts 


+ T 4S) (3-1) 


The block diagram form is shown in figure 3-5: 


am Geel 





Figure 3-5 The block diagram of an ideal 


PID controller 
The signals from the three branches in figure 3-5 may, at times, 
cancel each other partly. This situation can be visualized in 
figure 3-6 where the signs of the signals are considered. 
In the arrangement as shown in figure 3-5, whenever a signal 
in any of the three branches exceeds the output range of the op-amps, true 


representation of the wanted control action will be lost. If gain Ko is 
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Figure 3-6 Control action of different modes. 


positioned at the summing point instead of positioned before the 
three branches, the signals within the three branches will be 
considerably smaller (being only e(t)/T,, Ty e(t) and) escee(t)) .eewaito 
this arrangment, the output from the three branches will cancel each 
other partly, before being multiplied by the gain Ko. This arrangement 
will reduce the chances of saturation. 

However, the signal input to each branch must not be too small 
either, for noise considerations, The gain K. therefore is split 


into three parts i.e. 


K. = 0 x ies x ee (3-2) 
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where O<a<l, Ko! = 2, Koue = 10 

As mentioned in the section entitled "System description" in this 
chapter the gain K, is) lJimited to 20 in» this controller. With ao being 
represented by a potentiometer, Ko! and KO! areesetecoeesand lO 
respectively. Hence the gain K. can be varied by the potentiometer 


between values of zero and twenty. The format of the proportional 


gain is therefore decided and is shown in figure 3-7. 


E(S) 






we 
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Figure 3-7 The gain Ke is split into a, ie and KO 
to reduce the chan e of saturation 
The design of an active integrator is relatively easy compared 
toethe desiegn fot sardifferentiator jeeinathe wlatter, mtiltering of higher 
frequencies is absolutely necessary. A true differentiator can never 
be used because of its ever increasing gain for higher frequency 
components of noise and disturbances. Therefore, a close approximation 


rather than the ideal algorithm is all that can be achieved. 
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The circuit adopted is a passive lead network which, when 


combined with the branch of 5 wewLllabevvery close to the ideal 
a8 


aleorithm. sThis circuit is shown in) figure 3-8. 


Se 

Vo {e) 

Ry 

e(t) ee 
Re 
Figure 3-8 The passive lead network 
The transfer function of the circuit in figure 3-8 is: 
V(S) : SC[R, (R|+R,) + Ba Role eRe nas 
E(S) R)+R,+SC[R_ (RFR, ) + RR] 


The value of Rj and Ry» as later sections will show, are in the range 


of a few hundred ohms while R. is in the mega ohm range. Hence: 


V(5) . SC[R_ (R+R,)] + Bp 


E(S) ~SC[R_(R,+R,)] + Ry+R, 





=) 











S R + 1 
eRe oe: ey 
feCRedel | (a) 
ag R 
2 
CRe(RPaR *) R 
2 ale ee one im 2 , 
Let Ty = R and T, =c R. = qT; RR Then: 
D lee 
Vi() 3 Ts ce dk 
E(S) RR, 
Cres ) 
R 
2 
re T,s + 1 
it G_)) 
Ch po) aT 
d Ty 


Note that the ratio of ee over T; is a constant. The frequency 


response plot of equation (3-5), as given in figure 3-9, shows that 


= and — are the corner frequencies corresponding to the zero and 
d d 


pole location respectively. On the logarithmic abscissa, the two 
corner frequencies will be maintained at a fixed distance and hence 

a fixed gain difference, no matter in what position they appear along 
the abscissa. In this manner, whenever the position of the zero = 


is changed, the new position of the pole is also determined and any 


possibility of excessive gain at high frequencies is excluded. 


GAM ECEOG) 


G) (LOG) 





20 dB / DECADE 


Figure 3-9 The frequency response of equation (3-5) 
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The values of the lead network were chosen to be: 


9002 


Re 

R, = 1002 

Cae = at 

R. = 0.86 x M MegQ where 0 < M < 10 (3-6) 


Substituting equation (3-6) into (3-5) : 








Vi 6S) : ib BAe 1 
E(S) R)tR, 
Cleon) ) 
R 
2 
1+ 60 MS 


- “GQ ¥ 6 Ms) x 10 a. 
In equation (3-7) when Ta equals 60M seconds, Ta will be 6M seconds, 
in other words Ta = T,/10. As M varies between zero and ten,the derivative 
time constant Ty varies between zero and 600 seconds. With this set-— 
up, = will always remain one decade of frequency away from = A 
with a high frequency gain of 20 dB above the low frequency gain. 
As will be shown later, this distance of one decade and a gain 
different of 20 dB make the approximate algorithm close enough to the 


ideal one. The approximate algorithm in block diagram form is shown 


in figure 3-10. 
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E(S) 





Figure 3-10 The block diagram of the approximate 


PID algorithm 


The transfer function for the approximation is: 





1 -3res 
MOS) Aapsacia | d 1.5 
E(S) Cc fe sre ts TS 
Tk ge 12 os 
- d 1 
make eee sarase (3-8) 
1+795 


In most cases the time constant qT, used is of a rather small value 


[6leeethesvaltie of lf @which 1s only one tenth of 1_,) can be neglected. 


d d 
Therefore: 
M(S) _ il oe Abe l 
E(S) c Tere Tes T,§ 
ek eh Rae ee (3-9) 
sae d T,S 


It is clear from equation (3-9), that the algorithm adopted will give 
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a fairly independent parameter set as desired. To prove that the 
approximation is really a good one, the following analysis is carried 


Out.) Given the transfer function: 








fe) d di 
=a ( 1 ow) 
E(S) c T, 5+ 1 T,S 
where T,' = er (3-10) 
d Oma 
and the input e(t) is a ramp i.e. 
e(t) = At where A = constant 
A 
or E(S) = a (3-11) 
S 
Then: 
ae) Gy al 
ae ie Ne 
he) i, eres) aay a 
d ab S 
K ATS KA KA 
= Se aaa OF aaa a 3 (3-12) 
, Vo. S 
(Ty S+1)S (T Si.) o Ty 
Taking the inverse Laplace transform: 
ie eet G 
re Ti 
is d : i 
v(t) = K At +K AT (1 - e ) KAT , (l-e ) 
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(3-13) 
Equation (3-13) can be split into 3 parts: 
theseeaction 1s: Ye = K (At) 
: : ce At 
thewiy action. is: v, (t) =e Co 
: _ ioe 
and the Deaction is: Vy Ct) = K Tq: [O.9A(1 -e el )] 


Obviously, the P and I action in equation (3-13) are identical to 
the ideal P and I action. Only the D action in equation (3-13) shows 
the effect of approximation. The ideal D action for a ramp input 


would be: 
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(3-14) 
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To compare the ideal and the approximate D action, both of them are 


plotted in figure 3-11. 
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Figure 3-11 Comparison of the ideal and 


the approximate D action under a ramp input 
Figure 3-11 illustrates that at a time equal to T ,/10 the D action 
has reached 63% of its final value, and also the final value equals 
90% of the ideal value. 

Hence it becomes obvious that the algorithm adopted is a 
sufficiently close approximation. 

The integrator circuit is straightforward as shown in figure 


6-le.w either outputetrom the anteprator ci ecultsotetipure S-ireis: 


v(t) =a cx [2c dt 
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Figure 3-l2 The integrator circuit 


ax 


where v= 


ginal ©) < @ < Ib 
Al —_— _— 


Varying of integral time qT, is accomplished by means of a potentiometer. 
The op-amp chosen for integration must be selected with care to have 
fowl dritt and input oiiset) (16). 

The complete circuit for the PID controller proper is shown 
in figure 3-13. With regard to figure 3-13, the following points 
of interest should be mentioned. The value of M is designed to be 
adjustable in 20 steps. The first 10 steps vary M between 0.1 to ihe 
and the next 10 steps vary M from 1 to 10. The potentiometer at the 
summer is provided for padding the output to the range of 1 to 5 
volt. The zener diode acts as a voltage limiter. The switches 
around the integrator are there for the purpose of reset inhibition. 


Output current source and filtering 


The control voltage m'(t) generated by the PID summer will 


vie 
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PID controller proper 
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pass through a low pass filter, so that any frequency component above 
10 Hz will be attenuated. Mainly, the controller controls plants 
with large time constants in the range of minutes [6]. Therefore, 
frequencies above 10 Hz can be filtered out without affecting the 
relevant spectral characteristics of the control loop. At the same 
time noise and higher frequency components of disturbances are also 
reduced. 

The voltage to current conversion is achieved by the combin- 


ation of an op-amp and a medium power transistor shown in figure 
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Figure 3-14 The controller output driving circuit 


with low-pass filter 


The output of the op-amp directly drives the base of the power trans- 


istor. The true purpose is to keep V equal to Vn such that the 
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op-amp differential input voltage Neg will be zero. The value of 
Ro is chosen to be 1002. With 1 volt appearing across Ros there will 


be a current of 10 ma passing through it. When Beat aie Oli. 


E 
the maximum output, there will be 50 mA of current passing through 
Roe Since the current through Ro is I, (LAr )> approximately the 
Same amount of current will be forced through the collector load. 
Therefore, the intended purpose to convert output voltage 
Fane ingetlolert ecto t= VOLtsminto a 10 Marto +50 Masrange of current 
has been achieved. If a +40 volt supply is allowed while the output 


transistor is driven to saturation at 50 ma, the voltage value of 


Vemwiliebpes+soaVvoOlLeS. (understhe situation V 


RL =F+5)Vv0Olt); In other 


RE 
words, when the output current is at its maximum 50 ma, R can be 
chosen to have a value up to seven times the value of Ro before the 
transistor saturates. This is desirable, because there is no control 
over the transducer input impedance i.e. R- The best one can do is 
to set an upper limit for the transducer input impedance. In this 
case the +40 volt supply sets the upper limit at 7002. It should be 
remembered that this is a current sources therefore the transducer 
input impedance will not affect the current output. 

Note that the output from the summer is actually from -0.6 volt 
+ 6.5 volt instead of the specified +l to +5 volt range. This is 


dictated by the limiting diodes. When V = -0.6 volt, the output 


Ean 
of the op-amp becomes -15 volt in order to force enough reverse 


current through the reverse biased base-emitter junction. Understandably, 


this kind of operation will eventually damage the output circuit. 
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The remedy is to use a diode as a path by-passing the emitter- 
base junction when negative voltage appears across Ree This diode 
is also shown in figure 3-14. 

Reset Inhibiting 

The purpose of reset inhibiting was mentioned in Chapter II. 
A realization of it was also described. The method used was simply 
to monitor the output of the PID summer. Whenever the sum exceeds 
el thereousthembimi tS (touvolivoverle voltethe: integrator will be 
inhibited. 

First of all a direction detector is needed to indicate 
in which direction the controlled variable is moving. If it is 
rising towards maximum output, the PID summer output will be compared 
with +5 volt. Otherwise the summer output will be compared with +1 
volt. After the direction detector a set of special comparators 
is needed to check if either of the limits is exceeded. And finally, 
a set of solid state switches to carry out the inhibit or non-inhibit 
commands is required. 

For the direction detector, the circuit discussed in Appendix 
ietoechoseng 160] alne Inputestenal e, for it will be the negative 
error signal, -e(t), in the controller. Since -e(t) = c(t) - R(t), 
where R(t) is the fixed set point level, it is obvious that -e(t) can 
represent the controlled output c(t) for the purpose of direction 
detection. Note that the input e(t) is a signal that varies around 


the zero volt level. Hence the circuit in Appendix I will have eof 0° 
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The circuit of the direction detector is shown in figure 3-15. 
Whenever a major disturbance occurs, reset wind-up is 
possible. At the same time, the error signal will be large enough 
to cause a change of state in the direction detector. On the other 
hand, whenever the controlled variable approaches steady state, 
the error signal becomes very small. The threshold level in the 
circuit is designed so as to disallow the error signal under effectively 
steady state conditions to cause any change of state in the direction 
detector. It should be remembered that the integrator should never 
be inhibited under line-out (steady state) conditions. Moreover, 
noise of amplitudes smaller than the threshold level will not cause 
unwanted changes of state. 
The same level detector with hysteresis is used as the input 
stage of the signal comparators. Only this time Ore does not equal 


zero but equals either +5 or +1 volt. Suppose e equals +5 volts, 
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The asymmetrical values of e, are intended to act as the command 
signals for subsequent logic circuits. The logic circuits chosen 

are the RCA CD4016 COS-MOS switches. They are biased in such a manner 
that a +3 volt control signal is required. To obtain this +3 volt 


it was decided to pass the asymmetrical output e, through a gain 10 
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amplifier with a voltage limiter as mentioned in Appendix I. The 
complete comparator as well as the reset inhibiting circuit are 
shown in figure 3-15. The operation of the reset inhibiting circuit 
in figure 3-15 is roughly described in the following: Suppose there 
is a sudden positive step change in set point, and if it is assumed 
that the reset inhibiting is operating, there should be no overshoot 
in the controlled variable. The controlled signal, the error and 


set point changes are shown in figure 3-16. 
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Figure 3-16 The response of c(t) and e(t) 
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Theverror sienal, in) figure 3-16 is going positive. The positive error 
will cause point A and point B of figure 3-15 to become - 3 volt and 
+3 volt respectively. The values of voltage at points A and B will 
be maintained until an oppositely directed change of e(t) occurs. If 
the expected reset inhibiting action is going to take place, the 
error signal will not overshoot but will approach the new set point 
from one side only. Hence there will be no change of sign in e(t) 
and no change of state in the direction detector unless another 
disturbance should take place. At this very instant in time the 
voltage values at points A and B will cause switches 1 and 3 to close 
and switches 2 and 4 to open. The closing of switch 3 allows +5 
volts to appear as Soe at the positive end of the level comparator. 
At any moment when m(t), the control output, exceeds +5 volts, point 
c will become +3 volts and point D -3 volts. Recalling that the 
direction detector is still forcing switch 1 closed and switch 2 open, 
only +3 volt from point c will be allowed to pass through switch l. 
The result is that switch 5 in the integrator circuit will be closed 
and switch 6 will be open. Thus the integrator is inhibited. 

Once m(t) falls back below +5 volts, point c will reverse 
its sign to become -3 volts. Yet switch 1 is still closed. There- 
fore, -3 volts will pass through it and reactivate the integrator. 

When the error signal goes negative, switch 2 and 4 will be 


closed, making +l volt the e at the input of the level comparator. 


ref 


This time, only signals from point D can pass through switch 2 and 


control the integrator. Only when m(t) drops below +1 volt, will the 
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integrator be inhibited. 

There are three sets of mechanical switches involved with 
the integrator. They are: 

1) Reset inhibiting select switch. 

2) Reset on-off switch. 

3) Manual or automatic select switch. 


These switches together with the integrator are shown in figure 3-17. 
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When the reset inhibiting select switch is in the on 
position, the integrator is under the control of the reset inhibiting 
logic. When it is in the off position, a fixed -3 volts level or 
voltage appears and allows the integrator to function freely. 

There are occasions when reset action must be cancelled 
completely, such as in the case of pure PD operation. The rest on- 
off switch serves this purpose. 

The purpose of connecting the triple ganged auto manual select 
switch to the integrator is to provide bumpless transfer operation. 
When set in the manual position as shown in figure 3-17, the first 
of the triple ganged switches will provide +3 volt level to disable 
the integrator and at the same time turn the integrator into a gain 
l inverter. The second switch will provide -m,(t) ele vo Leto mehes an 
put of the integrator. Hence the positive manual control output minus 


one volt will appear at the output of the integrator. 


When switched back to automatic control, the integrator will 
be reactivated but with the initial value m,(t) —ievolt, stored an it. 
set point equals c(t) during the manual steady state,the automatic 
control output m(t) will be equal to the previous manual control 
out put m,(t) Thus no bump will occur during the transition. 

Note that all these mechanical switch outputs are tied to 
the same point. However, the +3 volts for the auto-manual switch 
and reset on-off switch are both derived from a voltage source, while 
the -3 volts for the reset inhibiting select switch comes from a 


voltage divider. Whenever +3 volt and -3 volt occur at the same time, 
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the +3 volt will always over-ride the -3 volts. In other words, 
whenever the reset action is at the off position (gives +3 volts) or 
when manual aes is selected (+3 volts also), the integrator 
ceases to be an integrator, regardless of the state of the other 
switches. 

Accessories 
A set of accessories on the front panel is required to 
select or indicate the proper functions. They are: 
1) Digital or manual set point select switch 
2) Manual set point knob 
3) Digital set-point switch set 
4) Automatic or manual control select switch 
5) Manual control knob 
6) Direct or indirect action select switch 
7) Reset inhibiting select switch 
8) Reset on-off switch 
9) Meter reading select switch, selects to show avesetepoine 
b) plant output, c) automatic control output, d) manual control 
output, e) zero check 
10) Digital set point read in switch 
11) Power on off switches for +15 volt and +40 volt supply 
12) Indicating lights for +15 and + 40 volt supply 
13) Current output connector: a female plug for cable connection 


at the back of the device 
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14) Voltage output connector: a female plug for monitoring purposes 
at the back of the device 
15) Input connector: a female plug for cable connection to receive 
feedback from the plant output 
16) Safety fuses: to limit power consumption from +15 volts and +40 volt 
supply 
17) P.I.D. parameter indicating dials 
The present circuit lay-out is purely for convenience reasons. 
The controller can be easily modified to a smaller size for industrial 
uses. The complete circuit, and the appearance of the device is shown 


in the following diagrams: 













uty Gataligiook 20% gulg ‘okie & ee ahi 
ibaa att ta asad 

avi vinit id. 'nabstanes ou hde ore qaly efiainas s -saseniies 2eee ae 
Subse” tots ana wre Appling? 4 - 7 

” o¢twy Blo geal “old geen thadees ‘yett 43 veded? -q70 108 TRE = 
ait ¢ | | a ee 
vba wettgs hed yeteeeteg 0,3, 7,032 a 


: Gs : “ae 
5 tyevews wot view eh, opieyel. aiede, peeesty 627.) : : 


: 
rlluse & so Sipe yp tees We ae rel lotaeee eee : 
, : st ; rq! wise SSS Ast ize Pee Ss ee |S) SJ >. Oa os » knew 


er 9a) ewes ofa pt 


2 az 
EXTERNAL SUPPLIES 









hy el ar ee | “ie 
meals Vet ge eet ce 
© 
POWER 
TRANSISTOR 
\NDIC ATION EeN4O\9 


LIGHT 






VOLTAGE 
REGULATOR 






INavwT2e8s 


t 
4 AMP FUSE = 





INDICATION 
LIGHT 


| 
| 
| STO 52 
| 
| 


SON, 


\NDOICATION 
LIGHT 





POWER 
TRANSISTOR 


2N4919 
VOLTAGE 


REGULATOR 
oe WE) 





FIG 3-18-O HOAs. “Sweiey ejli=veleiing 


| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 


| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
J 


cs ee 





vie Poe 


fr) p» 


nn 





(Tha iw 


ze FW ict 


iif 


els 
pts V 






READ IN " 
SWITCH RCA 
is) Is 
CD 4047 
+3 V 4 l2 
5 ONE 7 
SHOT 





GENERATOR !° 


FROM THE 
OTHER FOUR 
FLUP- FLOPS 





wuiosBrTes 432 


FAIRCHILD 
UA T22 


| 

|] DIGITAL SET | 
| POINT | 
& MORE SWITCHES ! 
| 

| 

| 





R ERE, REZRE, R E REREY 


4 MORE 
ELIPEELOeS 


FIG S3-\@b- DIGITAL SET POINT GENERATION. 





54 


"HSadOud BA NOYANOD AHL AO LINDABID AHL Del “D143 
Daa Nndwv = = 
SBASNVSAL Hoo! 
















W 4, SSaATIWwNA A\+ HS Law 
(i+™ Ww <j Wa a Ne 
( ) V Ss HO|\ = 
aool — 
Ae? 
NO 
ORIN YUN KO 4 Sng 

Hol DA SES HDLIMS 

es Cass 41DaNAaS 

ANAANEG yoo) = : <i WwalSw 

BOWLION fos fe QLOY acy gir WOl 
(yyy ob HOd = 
V Al- TOCYLNOD AAILVAISERG 
ee 2HescSe LINNO! 
Eee <8 Oot 
LNaxauND Se meicuc, 
ews i 
ama ols 
avon >> a Lee) 
< — 
] MO\ 
| 
ROT+ AAS TJOYULNOD AVUOBLNI = © AO} 


HOLVUHOBLNIA Jyico} YO Sol iL POSLNOS = = 


<G) 10d O01 rs 
Ol YOO 
4Oz7s <a) <a a. 
MOO} HOO\ Be G) ne _GBON 
a NO!LY} 
Gs) (5) AO! HOO! —— Sues a53uia! 
! 
CH) HOO! sno 


= NO\LOW LOAYIONI 


AOYULNOD (3) oe = 
<a) AN1IOd 18S $22 = 
ee ae WANUW HOS poles 

©) <a) Ol d = 

<3 HO! 

ITY OC »eo x 
OO! <<) 

oO! Hol ree 


SOO} mee \ vy, 
MZ 


ANGVD HOVEOS34 WO8S (QA)D LNdDLNO SSSDONd ——_+ 














COMPARATOR 





s)5) 





+ 
+ ev ine 
> << 
ee IooK «= 6CUBKK 
VK 
Fees IOK exe \OOK 5k 
‘ . iOK i 
(kK [Rca 3 
ICDAOIG 
_L SOK SWITCH 
= lz 
teen +IV 
DIRECTION DETECTOR 
Yeore T= at 
2 AN HSN ee ee 
CH TT — 
“1oOK — 
kK pies 
e A 
OOK 
Ueda: RESET 
INBIBITING 
e \OOK 
: en Oe 
ee een, NO te 002 
RESET UF 
COMPARATOR INHIBITING ii 5OK 
NY) I3V a 5 
ot<}- >t rt — 
i 5k = 


RES 


oO 
4+ 3 ACTION ON 


RESET 
AGTION © EE 


ElIG SA8- >) RESET INHIBITING Greet. 













eres a 
init 
| a i,’ —— aa 


aaa 


56 





—~\@ The Appearance of the near ideal P.I.D. controller 


Figure 3 





Chapter IV 


The Tuning of an Ideal Controller 


The existing tuning techniques are numerous. They are 
the product of different definitions of process and performance 
criteria [17, 18]. To name a few: there are Zeigler-Nichol's 
closed-loop and open-loop methods, Cohen-Coon's method, Harriot's 
method, three constraints method and integral of error methods. 
The co-existence of these methods proves two things: (1) control 
systems with nonlinear elements are difficult to understand 10 Te 
(2) optimal tuning brings marked improvement in process control, and 


hence the effort to search for it is justified [3]. 


A controller is not designed for the controlling of random 
noise, but rather for set point changes, disturbances and load changes 
[20]. The solving of a tuning problem usually starts from the proper 
characterization of the process itself and the choosing of a sound 
performance criterion. The object is to find a controller setting 
which will produce the desired performance for a particular plant. 

The characterization of a plant [19] 

The characterization of a plant falls into the category of 
open loop method. In contrast to open loop method, the closed loop 
method seeks to characterize the entire control loop. In practice 
it became obvious that the open loop is superior. Because the open 
loop method needs to upset the plant only once, while in the closed loop 
method the necessary information must be acquired through a number 


of upsetSof the entire control loop. In many cases the closed 
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loop method is not only cumbersome, but also intolerable. The 
open loop method requires information to be extracted from the process 
reaction curve. A process reaction curve can be produced in the 
following manner. 
1) Place the controller on manual operation and allow 
the system to reach steady state. 
2) With the control still in manual operation, impose a 
step change at the output of the controller. 
3) Record the response of the process output. 


A typical process reaction curve is shown in figure 4-l. 


MAGNITUDE OF 
PR@CESS OURmr oat 






SOS lan aa «ain oa Le 


TC 






POINT OF 


GO MAXIM UKMA 


ee i 


AO t 

a cateaia sae TIME (MINUTE) 
| 
| 


Figure 4-1 A process reaction curve 
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Different methods are available to approximate a process represented 
by the process reaction curve in figure 4-1. The most common approx- 
imation is a pure time delay plus a first order lag. One reason 

for the popularity of this approximation is that a real time delay 

of any duration can only be represented by a pure time delay. There 
is no other simple adequate approximation. A higher order lag plus 
pure time delay is also possible, but not necessary. The pure time 
delay plus first order lag approximation is easy to obtain and 
sufficiently accurate for most purposes. The approximated process 


transfer function is: 


mich 
ake 
Ae gost Get 
where oe = L. dead time delay 
ceded 180 alee ra &. time constant of a lst order 


system 

The parameters K, 8, and t can all be obtained from the process 
reaction curve as shown in figure 4-1. In figure 4-2 the reaction 
curve of a second order plant with dead time delay, and its approx- 
imation by a first order lag plus dead time delay are shown. 
Figure 4-2 reveals the closeness of the approximation to the true 
process reaction curve. 

or plants with an underdamped reaction curve, the first 
order lag plus dead time delay approximation is not suitable. Instead, 


a second order lag plus dead time delay model should be used. Note 
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PLANT OUTPUT 


ACTUAL PROCESS GS) = Se 
(Stal o.S SA) 







-{.91S 


APPROXIMATION G(S) “ip Te I 


Figure 4-2 The reaction curve of a 2nd order plant 


plus dead time delay and its approximation 


that this model can actually represent both underdamped and over- 
damped processes. This higher order approximation requires three 
Parameters to be extracted from the process reaction curve instead 
of two parameters as in the previous model. This extra bit of 
information insures a better approximation and as a result better 
controller tuning. 

The transfer function of the second order lag plus dead 


time delay model is: 


55 
G(S) = AS (4-2) 
Dorr ois $88 ee a 
2 W 
w n 
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There are a number of methods to extract 84? OF EE and K from the 
process reaction curve. Meyer's [21] method published in 1967 can 
be applied to both underdamp and overdamp systems with dead time 
delay. Meyer's method is good in the sense that the ease to apply 
is comparable to that of the first order model. 

The second order lag plus dead time method is a better approx- 
imation for higher order systems. However, in most cases the first 
order lag plus dead time delay approximation is sufficiently adequate. 
Unless a higher accuracy is needed or in the case of an underdamped 
reaction curve, the second order lag plus dead time model is seldom 
used. 

Choosing the performance criterion 

Most of the time the process under control is at steady 
State. Starting up, or set point changes are rare to encounter. 

Under the on-line control situation the most frequent happenings 
are load changes and disturbances. Hence the majority of control 
tunings are tuned for better performance under disturbance and load 
changes. The answer to what comprises the best performance, may 
well be a subjective oe. Most of the older tuning techniques were 
based on the so called quarter decay ratio as the performance criterion. 
The idea behind it is that, if the difference between each successive 
peak and the steady state value is one-fourth of the previous 
difference, the output will settle down very quickly. This situation 
is depicted in figure 4-3. 


However, the situation in figure 4-3 is only true for a pure 
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CONTROLLED vAaRisBtEe B _ 
Ge) = 





Figure 4-3 A response with a quarter decay characteristic 
second order system. In any other system, there is no guarantee 
that this decay ratio will hold true through all subsequent peaks. 
There is another problem with quarter decay ratio as a performance 
criterion. In a multi-parameter controller such as PI or PID 
controller, there are infinite combinations of parameters to produce 
a quarter decay response. It was based on this fact that arbitrary 
constraints were imposed on the original Zeigler-Nichol's open loop 
tuning method. It became known as the three constraints method bods 
The fallacy with quarter decay ratio is that, it relies upon a few 
points from the response curve to determine the goodness Olvascontreol. 
So is any other similar method bound to fail (e.g. frequency response 


method, also relies on a few points). 
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The ideal performance criterion should be a single figure 
that accounts for, not only a few points, but the entire closedloop 
response. Originally, the purpose of a feedback control system was 
to minimize the output error as a function of time after some upset 
has been introduced. Both the magnitude of error and the time 
over which the error is present, should contribute to the definition 
of the optimal control. Accordingly it is logical to propose a 
figure of merit to characterize the entire time response of the 


system as: 


co 


d = f reco, erae (4-3) 


fe) 
where F is a function of error and time. By integrating F with 
respect to time, a single number will be obtained, characteristic 
of the entire time response of the system. The criterion of 
performance: ® should have the property of being zero only if the 
error is zero at all times. This is clearly impossible after a 
disturbance ha$ been introduced. However, the smaller the value of 
@ is the better the performance of the system. Optimum performance 
can therefore be defined as that time response which gives the 
minimum value of @. Since the time response of a particular system 
is a function of the controller settings, the criterion of performance 
@ is also a function of the controller settings. The problem of 


minimizing © can be mathematically formulated as: 
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eles) (4-4) 
= 0 (4-5) 
= 0 (4-6) 
=a (4-7) 


where Roe T, 5 and T4 are the controller parameters for gain, integral 


time, and derivative time respectively. The simultaneous solution 


of equations (4-5) to (4-7) is the necessary condition to yield a 


controller setting which will cause an optimum performance. It is 


obvious then, 


that the use of a single figure of merit to characterize 


the system response results in a unique combination of parameters 


as the optimum setting. A number of possible definitions of 6 is 


listed below: 


wa 


[te dt = ITE integral of error x time (4-8) 
° 

aor 
je dt = ISE integral of error squared (4-9) 


00 
var dt = IAE integral of absolute error(4-10) 

a0 

t 


Je]dt = ITAE integral of absolute 
PTTAE vA 


error x time (4-11) 


00 
) = va te eit = ITSE integral of error square (4-12) 
ITSE f 


x time 
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Some of the above listed criteria are applied to a linear second order 


system: 


G(s) = +4 (4-13) 
Se Gt Pade ae Ae 


with a step function as the input variation. The resultant © value 


versus system parameter variation is plotted in figure 4-4[22]. 
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Figure 4-4 Integrals of error versus 
system parameter variation 
From figure 4-4, it is: clear that OT PAR? OAR and Pigg are the curves 


that have selectivity. Among the three, the performance criterion 
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OL TAR is the most selective one. When apply all different 4S to 


higher order systems, © still maintains its selective character 


ITAE 
while the other %'s wane into non-selective curves [22]. Hence, 

it can be concluded that OTAE is superior to the others, for both 

low and high order systems. In practice the most often used criteria 
ane: lon, stAReand sLrAr. 

The ISE criterion (integral of error squared) is relatively 
insensitive to small errors, while large errors contribute heavily 
to the value,of the integral. The use of ISE as a criterion will 
result in responses with small overshoots. However, long line-out 
time appears because small error occurring late in time do not 
contribute much to the csuuai= tenet. 

The IAE criterion (integral of absolute error) is more 
sensitive to smaller errors but less to large errors than in the case 
of ISE. Often the result of using this criterion is intermediate to 
thosevol LSivand  LTAR criteria. 

The ITAE criterion (integral of absolute error weighted by 
time) is insensitive to the initial unavoidable errors, but it penalizes 
heavily those error that occur late in time. Consequently, optimum 
response as defined by ITAE has the characteristic of short total 
response time and larger overshoots. 

Among these three criteria (ISE, IAE, ITAE), if sensitivity 
is not an issue, then it is left for the individual engineer to choose. 


The particularity of the system, or even personal taste can play a 


role in the choosing. 
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Tuning relationship based on integrals of error criteria 





With the process model properly defined and the performance 
criterion chosen, the focus should now be brought to the control 
loop itself. In figure 4-5, a simple control loop is shown with its 


controller blocks, process block and its inputs. 


ERROR ae DISTURBANCE 
SET POINT ENOL. 


CONTROLLED 
OUTPUT 














CONTROLLER CONTROL PLANT 


SIGNAL 






Figure 4-5 A simple feedback control loop 


The process block includes the dynamics of the measurement devices 
and that of the valves. Two input signals are shown, a change in 
setupoint, r(t) andea disturbance jsienal d(t). 

The transfer function for the control element is assumed to 


be represented by either: 


GS) = Ko P control (4-14) 
or G (Ss) =K (a+ a, PL control (4-15) 
(© Ts 
or CAG) ieee ae) ny Gekagal (4-16) 
c T,8 d 


Equations (4-14) to (4-16) represent the ideal controller algorithms, 
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with no lag in P or D modes and no interaction between all three 
modes. However, as was pointed out in Chapter III, these ideal 
algorithms can be very closely approximated by the controller proposed 
in this thesis. Hence, it can be claimed that any tuning relationship 
developed for an ideal controller, should also fit for this controller. 
This point will be shown to be true later in Chapter V. 

Assuming the process is a first order lag with dead time 
delay, and a PID controller is used, the response to a disturbance 


can be written as: 


-@s 

KeSte ~ 
2 eee ee Ls 
TS +Sc+ RK. e (T,S +S + i 





ce(S) = D(S) (4-17) 


EpeChesdvsturbance is@asunit@stepy then 


D(s) = (ChE 


and -6 8 
Ke 


c(S) - 
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(4-19) 
Re =) 


BE 


Making the transformation yw = 1S, the following equation results: 


_ ayn 
{6 
(4) Spi A SOR oe (4-20) 
te - Sur 
eae Ul ee eg (eae +ut>) 
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The above transformation is equivalent to introducing a dimensionless 


time t° defined by: 
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AS = 
ys 7 (4-21) 


Expressing equation (4-20) in time domain: 


c(t?) = K x(t’) (4-22) 
where 9% 
ee = 
=T ae 
K(t°) = [ aaa | 
fe) 
-—wH T 
2 i com 
uo tut KK e CG (WP anal i 
(4-23) 


Equation (4-23) indicates that the dynamic characteristic of the 


time response is a function of the dimensionless parameters: 


6 iE 6 
a KK. —aatia pa Of these parameters, only <= contains pure 


plant characteristics. The other parameters each contains a 
controller setting, and thus can be manipulated to obtain an optimal 
response. 


For a unit step disturbance, the error signal is given by: 
e(ty==— c(t) (4-24) 


Combining equations (4-24), (4-21), (4-22), (4-9) and using ISE as 


performance criteria, the result is: 


oo 
2 aye Oye a ee) 
oop = K of RN Es Ee (4-25) 


itll sis: ve. WOR) me lsaape 


; xe fl ig 4 










; ptiys itch lp, eS eieehtoay (bb) are ee 

itil of-Te peizsaw w @b ae ae 

(fo .<ahcanbtes paarta 30 “hie a «a oh 

PIN0> C945 23 ne 'oRns i] 5 ey } 5 sat ‘eun.sgpalipoaba 'yaktq: ~ 
_ vin lorie 60 gag apa ee pees satiexieee | 


RAS! 

N anmogagy 
® 

ban ie 7 


: " , hes 7 | 
‘yi povig oi Upaals actts ody (sone eee ite 2 0? _ i 7 
. a Arh _ ; 


= 
[ 
_ 


70 


The analytical calculation of the function eG2) involves the 
inverse transformation of equation (4-23). However, the presence of 
dead time in a control loop makes the analytical calculation of 

the inverse transformation extremely difficult. It seems reasonable 
then to consider numerical techniques to evaluate ISE in equation (4-25). 
In order to minimize the 6 in equation (4-25), the optimum 


ISE 


controller settings Ko T, and Ty have to be found. The multi- 
parameter optimization problem was solved by the use of a modified 
steepest ascent method known as optimum gradient. This optimization 
technique has been extensively discussed by Bekey [23]. 

Using integrals of error as performance criteria A.M. Lopez 
set up a control loop simulation and a three dimensional a ORT EE 
program on an IBM 7040 digital computer. The result was a set of 

tuning relationship graphs and equations. Those tuning graphs and 
equations that are related to this thesis are listed in Appendix III. 

With the understanding that the designed controller should 
be a very close approximation to the ideal algorithm, A.M. Lopez's 
tuning graphs and equations should also be directly applicable to 
the controller designed. These graphs and equations will surely enable 
one to determine the optimum settings rapidly. They also eliminate 
the numerous trials necessary in the older tunning techniques. 


For more information on tuning please see references 24, 25 
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Chapter V 


Results and Conclusions 





In order to test the designed controller, the complete 
control loop must be set up. The plant itself with dead time delay 
were simulated by a hybrid system, consisting of a PDP8 digital | 
computer, a EAI TR-48 analog computer and a REDCOR 610 linkage system. 


A detailed diagram is shown in figure 5-l. 


ia eae ors ras ra age ay mye) ‘aie | 
| 0(S) : 
| | REDCOR le 
IR& E(S) M(S) SIO D 
|— MGS es) : A/D DIGITAL. 
| Ws | mys COMPUTER 
| — Te | : CONVERTER 
| l 
| | 
| | 

| 
[Sa alies | 86 Per © MG OLEERE amy 

EAT TR-ABSB 
ANALOG 
(Sits) COMPUTER 
c(s) 


| HYBRID Ne 


SYSTEM PLANT SIMULATION 


Ne et cca, ally Aa Sa Tle egetint peace ae Ea 


Figure 5-1 The complete control loop with the plant and dead time 


simulated by a hybrid system 
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From figure 5-1 it is clear that the D/A and A/D converter together 
with the digital computer functions as a dead time generator, while 
the analog computer simulates the plant proper. The control signal 


plus disturbance: 


M (8) = M(S) + D(S) (6=1) 


is sampled by anA/D converter at a certain rate. This rate is 
BO72 


fe) 
have been given to the capability of the hybrid system and the 





determined to be times/sec (Appendix A II). Considerations 
desirability of higher rate of sampling. 

The sampled analog signal is converted into binary code and 
stored in a specified memory area of the digital computer. After an 
units of time, this binary code will be brought out and converted 
back to an analog signal. The reconverted analog signal will then 
be sent to the simulated plant in the analog computer. The detailed 
explanation, flow chart and program are listed in Appendix A II. 

The following plant models were used during the test of the 


designed controller: 
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For all practical purposes, these two models can represent most of 
the processes. Plants without dead time delay are not considered here 
in order to shorten this presentation. Furthermore, only PID controlling 
form will be used. 
The original objective of the proposed controller is to 
build a near ideal PID controller, so that A.M. Lopez's optimum 
tuning chart can be utilized. Whether this design objective is achieved 
Ob noes canebe on by comparing the results obtained from different 


tunings. In all subsequent tests, the IAE performance criterion 


is used. 
First order plant with dead time delay a 
fo) 
| SA rs : 
The chosen plant parameters for G,(S) a ace alan ot me 
K=1 
and 
t = 2.1 min = 126 sec 
IS) 
Three arbitrary = values were chosen to test the designed controller. 
8 8 
They are : ae On 0 SbandeOe Oo mel stinge uw sel 7Ousce, = 0.2 means 
aC ) 
6 = 25,2egec, —--=-0.5-means 6 = 63 sec, and — =-0.8 means @ -= 100.8 
O c O i Oo 


8 
sec. For each of the = values, Figure AII-I contains the optimum 


P, I, D settings. They shall be designated as Ko? Tee and Tos 
The controller will then be tuned according to these optimum values. 
Set point r(t) will be held constant at 2.5 volts, a mid point value 


of the output span. Then the control loop will be allowed to come 


to steady state. When steady state arrives, a disturbance of 1 volt 


e008 


wiian 0,0 <= Bs tal Bee ie RO. shyt os wee 
é oo = Oj see £3, & eau ae 
- muni igo gia acbprincn MbEn veg ee bt 


4 
4 


















Jeon 3 dag. Ges afeden vk tide 
ob hal a ta wv Poth ee elise tome 
AF ¥ing isisue han Sebice ep alas et tnd .a? 7 
.beeu ot Liaw ; 
of af galiorgoe beeweatdq SAS eyvbipel de Jani gere are . ; 
i aon Mah Neel aie ipaLbeuance Ult Losbe spam, eBid ’ 
wertaauril yom: waa ‘aoe ‘besdbeso «sd 17> reds gulag 
iti hbetletde: etteaee nb g0izequian wd gets wd ean gon'ae - 
nn MA et ede Pnqiges tye ofa af gue a 

oa poe st 
yals> outs dash fate smoke 3 rae ab, a 


BS et i oe wats io €IStcts tag 36834 wegetis oct 
a 
4 m2 a A is 
oP 
| ee 


ae 5a Coit > ee 35 * FP: 
a : ‘ & 
aru b Aa re js2 O23 \TmiOHe gaa smnQ ei, 


SRT PS na tah oY wate _ gots. ~ 
etic aanial oe oe a 


a: = 


74 


magnitude (a substantial one) will be imposed onto the control out- 
DuCs me lhe control loop is upsetted and the output response will be 
recorded. This response should be the optimum one. 

Then each of the three parameters will take turns to vary 
above and below the optimum values while the other two parameters 
remain at their respective optimum values. Therefore, a total of 
six variations will be obtained. Each of these variations should 
bring an inferior response to the supposedly optimum response. 

The test was carried out and the results are listed in 


Figure 5-2 to 5-10. 





Second order plant with dead time delay 9 S 
fe) 
In the transfer function G,(S) = ee »K and 
2 S Des 
——~ + + 1 

2 Ww) 

W n 
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wo are chosen to be: 


and 


Two arbitrary values of 8° n were chosen. They are: 88h = 1.1 and 
8 8A = 4.0. The damping ratio § is allowed to be either 1 or 0.5 to 
represent overdamp or underdamp systems. Therefore, four different 
sets of plant characteristics are available. Optimum PID settings for 
them can be obtained from figure AIII - 2 to AIII - 10. 

The controller was tuned according to the optimum settings 


and its deviations. The supposedly optimum and non-optimum responses 


were recorded in figure 5-11 to 5-22. 
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Test of reset inhibiting action 


Extensive trial of reset inhibition was carried out. In 
figure 5-23 the responses of a particular first order plant with 
dead time under a set point change are shown. Between the two curves, 


one is under the effect of reset inhibition. 


It is obvious from figure 5-23 that the designed controller 


has an excellent reset inhibiting feature. 


; : F ; i 
Extreme situations were also tried. With KK. = 2, 7 = 10 
al 
repeats /minute and T, = 0 minutes, the same first order plant in 


d 


figure 5-23 was put under a set point change and reset inhibition. 
Overshoot was found in the response. This was because of the extremely 
: i al es 
high value of TO Even when error started to diminish, the 
i 
integrator still charged up fast enough to cause an overshoot. 

However, an extreme parameter combination like this case is 
unlikely to happen in real situations. The purpose of integral action 
is to provide offset free on-line control. Its capacity to cancel 
offset is determined by the integrator output range instead of the 

i : : 
integral time: T,- The parameter T) only determines the integrator's 
Z 
speed to acquire its offset cancelling capacity. With critical processes 
such as batch sintering, batch brazing, diffusion, crystal growing 
etc., manual start up (very slow) is used to avoid overshoot. This 


suggests the speed to arrive at steady state is of lower priority 


than to avoid overshoot. 
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Therefore, when reset inhibition is needed, it is recommended 


to have the parameter _ set towards lower values. The relationship 
ah 

between Ko Ty and Tt. for an optimum control under reset inhibition, 
is complicated enough to warrant another investigation. 
Test of bumpless transfer 

Various processes were allowed to come to steady state 
under manual control. With the set point value set equal to output 
value c(t), the control was suddenly returned to automatic. The 
resultant responses observed were bumpless. 
Conclusion 

The objective to build a near ideal PID controller was 
achieved with negligible interrelationship between the three parameters. 
The PID algorithm adopted in the design was exactly the same algorithm 
which A.M. Lopez used for his optimum tuning study. Hence, Lopez's 
optimum tuning graphs can be readily applicable to this particular 
PID controller. As the three parameters are practically independent 
of each other, tuning the controller becomes a straightforward task. 
The recordings in figure 5-2 to 5-22 proved that this controller, 
when used in conjunction with the optimum tuning graphs, produces 
superior responses than those supposedly non optimally tuned results. 

The low drift integrated circuits used in this design, make 
it possible to abandon the DC to AC and AC to DC processing of signal 
as was necessary in the old days. 


The reset inhibiting feature is very satisfactory except in 


extreme cases. But as explained in the previous section, most of the 
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time these extreme cases do not exist. A lower Tt value can 
ak 


always circumvent this problem without impairing the controller 
capability. Themrecee inhibiting feature is expected to increase 
this controller's versatility. 

The manual control functions properly. The transfer from 
manual to automatic control is always bumpless. Other features 
such as selected meter read-out, analog and digital set point adjust- 
ment, power supply indicators, etc. are all meant to raetease the 


convenience in using this unit. 
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Appendix I 


Level detector with Hysterisis 


The circuit of the level detector with hysteresisis shown 


in figure AI-l. 





Figure AI-1 Level detector with hysteresis 


The output voltage, 


= +t + 
A VE Pee 
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Substitute equation 1 into equation 2 


sj R, Ss Ry 3 ref 
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(1) 


(2) 







1 ibe. 


ele vope\ 8 aesoused teva’! 


7 a 


ed 


‘ | herrea dl‘. ete lates eft 70 ston oft 
 «t-Ta soon ad 





ee 
ne 
— 

J 

a 


jew 3 akon oath helen 1 aa i = 


104 


+e (3) 


Substitute equation 3 into equation l 


Se 
et —_ 
te ‘= Ve = wy Ry ue fref 
Or 
R, +R 
4 3 
Ge 
ey fe Vv, ( R, ) + Cree (4) 


finitely e. is smaller than ene the current i, will 


flow from ey towards Be and e.- Vy under this situation will take 


on a positive sign. Using equation 3: 


s 
= + — + 
“54 a RD “ref 
and 
i ae OR 
Sivan ie XC 2 2 +e (5) 
fo) L Ry ref 
a3 
Tfeemerowseandwexceedssthesvalucseuwesat sate » the current 
a Sj ih Ry, ref 
i, will) reverse dts direction and Vy. will take on a negative sign. In 


other words, 





R 
x 3 
“sj a Ry, uF “ref 
R_+R 
= oie nah 
Soe Me ( R hat ref (6) 
4 
From this moment on, so long as ey is greater than the value De 
R 
(re =- vi. + Cef)? the output will remain at the value of 
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R,+R 
ame 
¥ a Ry et “ref’ 
LE e, varies its direction and becomes less than aH 
R 
3 
a et = A F é 
( Ve R, eof ome the output ey will switch back to the value: 
RoR R 
+ . i Poca 
vi, ( R i) Se e ef So long as e, is smaller than + Vi. R “f at 


4 4 
the output value will remain the same. The complete situation is 


depicted in figure AI-2. 
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Figure AI-2 Voltage output versus voltage input of the 
level detector circuit with hysterisis 
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level. The value e. has to exceed cnee plus or minus this threshold 


Figure AI-2 shows that the value |V 
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changed, the output will be kept in that state, unless e, reverses 
its direction and exceeds the opposite threshold point. With this 
hysterisis, the changes of states at the output will be the results 
of substantial changes of e.. Noise and errors of values smaller 
than the threshold will not cause unwanted changes of state. 

One special case of this circuit occurs when the reference 


voltage e. 


oe equals zero. The circuit becomes a zero crossing 


detector with hysterisis. Its voltage output and input is shown in 


figure AI-3. 
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Figure AI-3 Cr of the zero crossing detector 


with hysterisis 


Note the symmetry in figure AI-3. The two states of output e, are 


equal in magnitude but opposite in sign. This symmetrical output ae 
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can be used as logic commands on symmetrical COS-MOS switches such 
as RCA's CD4016. 


In general, when e does not equal zero, unsymmetrical 


ref 
output will occur in figure AI-2. If a symmetrical output is needed 
one can use a special circuit to transform the output in figure 


AI-2 into a symmetrical one. This transforming circuit is simply 


a voltage amplifier with voltage limiters. It is shown in figure 
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Figure AI-4 Transforming an unsymmetrical signal into 


a symmetrical one 
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Appendix II 
The Generation of Dead Time 
If a signal x(t) passes through a device which is nothing 
but a pure time delay generator, then the output y(t) should be 
exactly the same as input x(t), only delayed by tbe units of time. 


This situation is depicted in figure AII-l. 


MAGNITUDE 







(Et) INPUT 


y(t) OUTPUT 


Figure ALI-1 y(t) = x(t-6), T= sampling period 


9 = 3072 x T 
fe) 
This pure time delay generator can be simulated as in 


figure AII-2: ; 
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Figure AII-2 A simulated pure time delay generator 
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The control of the dead time delay generator is done by a computer 
program. At the starting point, the control program will clear all 
the core content in a specified memory area. Then, starting from 
the first address in that memory area, the program brings out the 
content of that address. The content is converted to an analog 
signal and sent out as output y(t). Since the contents of the 
memory is cleared initially, y(t) is zero at this instant. Next, 

the program will command the A/D converter to bring in the first 
sample of input x(t) in digital format. This first sample of input 
is subsequently stored into the first address of the memory area. 
After this, the program will wait for '"T" seconds, before it repeats 
the above sequence for address number two. The total number of 
addresses in the specified memory area is 3072. When the program 
comes to the last address (the 3072th address), T x 3072 seconds have 
elapsed. Then the program will go back to the first address and send 
out its content. However, the content in the first address is actually 
the input signal sampled at T x 3072 seconds ago. The program will 
run through the memory area repeatedly. If am UsmspDecd ated amtne 
value of "T" can be programed in such a way that T x 3072 equals oe 
units of time. In other words, for any ol value, one should sample it 
into 3072 intervals. The number 3072 is three quarters of the total 
memory address in a 4K PDP8 computer. Also, it is an extremely high 
rate of sampling for this particular purpose where oe usually 

exceeds tens of seconds. The flow and format of information is shown 
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SAMPLING 
A/D 
p= DIA = 

ANALOG DEVICE CONVERTER FITTER ANALOG 
SIGNAL SIGNAL 
INPUT OUTPUT 

action action 

SOre, Bove 


Note: The primed action comes 
right behind the unprimed 
action. Also actions of 
consecutive number are spaced 
"T" seconds apart. 


Figure AII-2 The flow and format of information in 


a pure time delay generator 
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THE DEAD TIME SIMULATION PROGRAM 


ATOD 8231 
CLEAR 6266 
CONST 6267 
COUNT 64278 
DACHN 6272 
DTOA G228 
IMDEX 6611 
INDEX 6818 
KONST 6273 
KONST1 @274 
KOUNT 96275 
KOUNT 1 9276 
EOPI 258 
NINDE 6266 
NULO G21 
OLOP S217 
TABLE 1888 
TEST 9269 
TTABLE 6271 


G16 BG0S 
G11 GGG 
C266 7388 
261 1267 
202 32786 
203 1271 
G2eg4 3616 
G2G5 3416 
$296 2276 
267 S265 
Geng” Mews 
@211 3916 
Yee shede 
8213 3811 
214 7386 
$215 T4G4 
216 3274 
217 7360 
G2eeG 1416 
221 6462 
G222 6464 
§223 7208 
fey -e a Wo fr 
0225 6451 
G226 6454 
G2e7 6441 


0236 S2e7 


*1G 

INDEXs @ 

IMDEXs, @ 

*290 
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Appendix III 
OPTIMAL TUNING GRAPHS BASED ON INTEGRALS OF ERROR CRITERIA 


ARO, Lopez in his Ph.D. dissertation had set up a multi- 
parameter optimization program to find out the optimum settings of 
a controller for particular control loops under specified performance 
Crire On. 

The performance criteria chosen were of the integral of 
error type. The plants chosen could be either first order plant with 
dead time delay or second order plant with dead time delay. The 
multiparameter optimization technique used was a modified steepest 
ascent method known as optimum gradient Peo 

Both the control loop simulation and the optimization program 
were carried out on an IBM 7040 digital computer. The result was 
a set of tuning relationship graphs and equations. 


For a first order plant with dead time delay: 


G(S) = 


the optimum PID controller settings for ISE, IAE and ITAE performance 
criteria are shown in figure AIII-1. Im addition, equations in table 
AIII-1 were developed to fit the curves in figure AIII-l. Fora 


second order plant with dead time delay: 


~ 655 
Ke 
2) 
pamee Sane 
p. Ww 
w) n 
n 


G(S) = 


the optimum PID controller settings for ISE, IAE and ITAE performance 


criteria are shown in figure AIII-2 to figure ALLE LO 
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For the optimum controller settings for one mode and two 


mode controlling and the computer programs please refer to A.M. Lopez's 


Ph.D. dissertation [6]. 


All graphs and tables in this Appendix are extracted for A.M. Lopez's Ph.D. 


thesis. 











PID Controller Settings to Minimize Error Integrals, 


Figure AMl-{ 
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TABLE AIL-| 


Tuning Equations for PID Controllers. 


Controller Algorithm: 


Tuning Equations: 


Criterion . a 


ISE 1.495 
IAE 12535 
ITAE 1.357 


-0.945 


-0.921 


-0.947 


m(t) 


x, (i+ oe + Typ )e(t) 
- «(2) 
(2) 
Cay 


d e £ 
-0.771 0.560 1.006 
ee ee cy oe 
-0.738 0.381 0.995 
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Figure All-2 Optimum Gain for PID Control. ISE Criterion, 
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Figure AJJ-3 Optimum Reset Setting for PID Control. 
ISE Criterion, 
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Optimum Gain for PID Control. IAE Criterion. 


| Figure Ail-5 
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ITAE Criterion. 


Optimum Gain for PID Control, 


Figure A w- a 
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